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Abstract — The encoding complexity of a general (en,ek) quasi- 
cyclic code is 0{e^{n — k)k). This paper presents a novel low- 
complexity encoding algorithm for quasi-cyclic (QC) codes based 
on matrix transformation. First, a message vector is encoded 
into a transformed codeword in the transform domain. Then, 
the transmitted codeword is obtained from the transformed 
codeword by the inverse Galois Fourier transform. For binary 
QC codes, a simple and fast mapping is required to post-process 
the transformed codeword such that the transmitted codeword 
is binary as well. The complexity of our proposed encoding 
algorithm is 0{e{n—k)k) symbol operations for non-binary codes 
and 0(e(n — fc)fclog2e) bit operations for binary codes. These 
complexities are much lower than their traditional counterpart 
0{e {n~k)k). For example, our complexity of encoding a 64-ary 
(4095,2160) QC code is only 1.59% of that of traditional encoding, 
and our complexities of encoding the binary (4095, 2160) and 
(8176, 7154) QC codes are respectively 9.52% and 1.77% of those 
of traditional encoding. We also study the application of our 
low-complexity encoding algorithm to one of the most important 
subclasses of QC codes, namely QC-LDPC codes, especially when 
their parity-check matrices are rank deficient. 

Index Terms — Quasi-cyclic codes, LDPC codes, encoding com- 
plexity, redundant rows, matrix transformation, and the Galois 
Fourier transform. 



I. Introduction 

Quasi-cyclic (QC) codes (T\ are an important class of 
linear error-correcting codes in both coding theory and their 
applications. These codes can asymptotically approach the 
Varshamov-Gilbert bound f2\. Moreover, their partial cyclic 
structure simplifies their encoding and decoding implementa- 
tions by using simple shift registers and logic circuits ||3J. In 
recent years, research on QC codes has focused on one of their 
subclasses, known as QC low-density parity-check (LDPC) 
codes Bl- lfm . which have been shown to perform as well as 
other types of LDPC codes in most applications. QC-LDPC 
codes have advantages over other types of LDPC codes in 
hardware implementation of encoding [12] and decoding [5], 
|[rr|. Thus, most LDPC codes adopted as standard codes for 
various next-generation communication and storage systems 
are QC. 

One important development of QC codes is the introduction 
of matrix transformation via the Galois Fourier transform ifTSl . 

This work was supported by National Natural Science Foundation of China 
(61201156). 



lfT4l . In the Fourier transform domain, an array of circulants 
is specified by a diagonal matrix over a finite field. It has 
been shown to be amenable to analysis and construction of 
some QC-LDPC codes flSl . Furthermore, the authors of |1T6) 
expanded the analyses of ranks and row-redundancies into 
most QC-LDPC codes. 

QC codes generally are encoded by multiplying a message 
vector m of length ek with an ek x en generator matrix G, 

where G usually is systematic f3| with G = [I : P] and 
I being the identity matrix — in the rest of the paper, G is 
systematic if not stated otherwise. There are two issues with 
the implementation of QC code systems. First, the generator 
matrix G is usually not sparse, so it requires a large number, 
i.e., e{n — k)k, of memory units to store P. Second, although 
encoding of QC codes can be partially parallelized so that 
the computation units are reduced by a factor of e, the total 
number of symbol operations is still e^{n — k)k, which is the 
same as that for general linear codes. 

In this paper, we propose to encode QC codes in the Fourier 
domain rather than using direct multiplications in the symbol 
domain. We are motivated by the fact that the row space of 
the transformed generator matrix is the null space of the trans- 
formed parity-check matrix. Moreover, the ek x en transformed 
generator matrix is an e x e diagonal array of fc x n matrices. 
Consequently, encoding in the transform domain (ETD) is 
achieved by e times encoding of message vectors of length k 
with k X n generator matrices (rather than encoding of an ek 
message vector with an ek x en generator matrix). We call the 
resulting vector after ETD the transformed codeword. The final 
transmitted QC codeword is obtained by permutations and n 
times inverse Galois Fourier transforms from the transformed 
codeword. Thus, the computational complexity of ETD is 
0{e{n — k)k), which is much lower than that of traditional 
encoding. The memory consumption of ETD is 0{e{n — k)k) 
symbols, which is the same as that of traditional encoding. 
For example, we show that the computational complexity of 
ETD of a 64-ary (4095,2160) QC code is only 1.59% of that 
of traditional encoding. 

Since binary QC codes are used in many applications, we 
carry out a detailed study of ETD of binary QC codes. For any 
binary QC code, its transformed generator matrix satisfies the 
conjugacy constraint ifTTl . ifTSl . but the transmitted codeword 
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from ETD usually is not binary. Thus, it costs several bits 
to transmit a code symbol, resulting in lower code rate. To 
make the transmitted codeword binary, we devise a simple 
and fast mapping to post-process the transformed codeword so 
that the new transformed codeword still satisfies the conjugacy 
constraint. The post-processing step consisting of a mapping 
with bases of subfields is the key in the binary case. The other 
steps of ETD of binary QC codes are the same as those for 
the non-binary case. Furthermore, if we take advantages of 
the conjugacy constraint on the transformed generator matrix, 
the computational complexity of ETD can be reduced to 
0(e(n — fc)fc log2 e). Its memory consumption is 0{e{n — k)k), 
the same as in traditional one. We show that the computation 
complexities of ETD of the binary (4095, 2160) and (8176, 
7154) QC codes are respectively 9.52% and 1.77% of those 
of traditional encoding. 

ETD readily applies to QC-LDPC codes with full-rank 
parity-check matrices. However, there are many QC-LDPC 
codes, especially algebraic LDPC codes Q, IE], IHl, ifTSl - 
||23l, whose parity-check matrices are abundant in redundant 
rows, i.e., rank deficient. By carefully constructing the trans- 
formed generator matrix according to the rank of the diagonal 
matrices on the transformed parity-check matrix, we show that 
our proposed ETD algorithm works for such QC-LDPC (or 
QC) codes as well. 

The rest of this paper is organized as follows. Section 
II introduces QC codes and the matrix transformation. In 
Section III, we present encoding of QC codes in the transform 
domain and study its computational complexity and memory 
consumption. Section IV focuses on ETD of binary QC codes 
and its simplification due to the conjugacy constraint. The key 
step in the binary case that involves construction of a post- 
processing mapping using subfield bases is derived. Section V 
is concerned with QC-LDPC codes, especially whose parity- 
check matrices are rank deficient. Section VI concludes the 
paper. 



II. Introduction of Matrix Transformation 

In this section, we briefly describe QC codes and matrix 
transformations of their generator matrices and parity-check 
matrices and refer readers to lTl4l . HI 51 for more details. 

Let W ~ [wij], < z,j < e, be an e X e circulant matrix 
over GF(g), i.e., if every row is a cyclic-shift (one place to 
right) of the row above it, including end-around. Then, we 
write W — circ{w), where w is its top row, called the 
generator of W ll24l . A binary (en,ek) QC code of length 
en is given by a fc x n array G of e x e circulant matrices (or 
simply circulants), named the generator matrix. 



G 



Wo.o 
Wi.o 



Wo,i 
Wi,i 



w 



fc-1.0 



w 



fc-1,1 



Wo.„- 



(1) 



or the null space of an {n — k) x n array H — [A^ j], parity- 
check matrix, of circulants of the same size, where G H^ = 
and Ai.j = circ{Biij) is a circulant of length e. 



Consider a finite field of 2"" elements GF(2''), i.e., q = 2*". 
Let a be an element of GF(2'') with order e and w be an 
e-tuple over GF(2^), where e is a factor of 2'' — 1. Here 
we consider only the case e = 2'' — 1. The other cases are 
similar The Galois Fourier transform of the e-tuple w fTTl 
over GF(2''), denoted by J^[w], is given by the e-tuple d over 
GF(2'') whose t-th component, dt, for < t < e, is given 
by dt = wq + wia^* + W2a~^* + • • • + Wg-ia^'-'^""'"''*- The 
vector w, which is called the inverse Fourier transform of the 
vector d, denoted by w = J^^^[d], can be reconstructed from 
do + dia' + d2a^' H h dg-ia^^-^)', for < / < e. 



Wl 



Consider an e x e circulant W = circ{w). Define two e x e 
matrices over GF(2'') as follows: V = [a^'-'] and V^^ = 
[a*-'], < i,j < e. Both matrices, V and V^^, known as 
Vandermonde matrices ifTTl ll25l . are non-singular. Moreover, 
V^^ is the inverse of V and vice versa. Taking the matrix 
product V^^WV, we obtain the following e x e diagonal 
matrix over GF(2''): 



W-^ ^ V-^WV = diag{do,di, 



(2) 



where the diagonal vector (do, di, rfe-i) is the Fourier 
transform of the generator w of W. The diagonal matrix 
W-^ — V^^ WV is referred to as the Fourier transform of the 
circulant W. If w is an e-tuple over GF(2), the components 
must satisfy the following constraint Iil71 



(3) 



for < i < e, where {2t)e denotes the nonnegative integer 
less than e and is congruent to 2t modulo e. This condition is 
known as the conjugacy constraint, which is the key constraint 
of the binary case in Section IV. 

Let k and n be two positive integers. Let G = [W^j], 
< i < fc, < J < n, be an fc X n array of e x e circulants 
Wi j over GF(2''), where 
W 

and n~^{n) as an n x n diagonal array of V^^ 
n{k) = dmg(V, ...,V), 

k 

n-^{n) = diag{Y-\...,V-^ 



J.J is the generator of the circulant 
ij. Next, we define il{k) as a fc x fc diagonal array of V's 

s, 



(4) 



Then the Fourier transform of G 

n-^ik)Gn{n) = [Wf^.], where 
X e diagonal matrix 



IS an 



IS given 
with diagonal 



as 
V 



G^ 



(5) 



vector 



z,j,0 7 1 1 



di j e-i), which is the Fourier transform of 



the generator j of j. 

Define the following index sequences: for < i,i < e, 

T'row,i = [i,e + i,...,{k - l)e + i] and ttcoI^j = [j,e + 

J, (tZ l)e -\- j^. Let IT row — [^ro-ujjO; '^row^l ; ■ ■ • ? '^rou?,e— l] 

and TTcoi = [7rco;,o,7rco/a,---,7rco;,e-i]- Then iVrow gives a 
permutation of the indices of the rows of G"'^ while tTcoI 
represents a permutation of the columns of G-^. Their reverse 
permutations are denoted by ■Krow ^^'^ ^coi^ respectively. We 
define the permutation tt that performs both row and column 
permutations. Its reverse permutation is denoted by tt^^. By 
the permutation tt, G-^ results in the following e x e diagonal 
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array of k x n matrices over GF(2''), 
G^'^ = dzag(Do,Di, 



(6) 



The transformation from G to G-^ '^ through G-^ is reversible. 
The reverse process is called the inverse matrix transformation, 
denoted by {T^^,tt^^}. 

If the array G of circulants and zero matrices (ZM) is over 
GF(2) C GF(2''), the matrices on the main diagonal of the 
array G-^'^ satisfy the conjugacy constraint ifTTl . 



D 



(2t)e 



D 



o2 



(7) 



i.e., the entry at location of D(2t)^ is the square of the 
entry at location of Df. We call the matrix D(2t)^ a 

conjugate matrix of Df . Following the definition of conjugate 
matrix, we can group all the matrices on the main diagonal 
Di, z = 0, 1, . . . , e — 1, into conjugacy classes. Let A be the 
number of distinct conjugacy classes and ^fo, 'Si, . . . , 
represent these classes, with 



,D 



o2 



,D 



02"^ 



U 

(8) 

where r)i is the number of matrices in the conjugacy class 
^i, i.e., rji is the smallest nonnegative integer such that 
{2^Hi)e — ti, and ti is the smallest number in the subscripts 
of the conjugate matrices in ^fj. The member matrix Df. is 
called the representative of the conjugacy class 'Si. 

III. Encoding of QC Codes in the Transformed 
Domain 

In this section, we present our proposed low-complexity 
ETD algorithm of QC codes. First, we derive ETD based 
on orthogonality fT5\ of the transformed parity-check matrix 
and the transformed generator matrix. Then, we compare the 
computational complexity and memory consumption of ETD 
with those of traditional encoding. 

Consider an (en, efc) QC code C over GF(2'') defined by the 
ek X en generator matrix G in ([T]), which consists of circulants 
of size e. Suppose that a is a primitive element in GF(2''), 
i.e., e = g — 1 = 2'' — 1. By the matrix transformation (|6), the 
generator matrix results in the transformed generator matrix 
Qj^,7r Qygj. Qp(2'"). It can be employed to encode a message 
vector m of ek bits into a transformed codeword c-^''^ of 
length en bits. Since 

G • H"^ n{k) ■ n-\k) ■ G ■ n{n) ■ n-\n) • h"^ = o, 

where H is the parity-check matrix of C, then 

n{k) -G^ ■n-\n) -11^ = 0. (9) 

Because il{k) is a nonsingular matrix, 

G^ -n-^n) -11^ ^0. (10) 

Furthermore, permuting the rows of G^^ by iTrow results in 
G-^.^™-, then 



1 J-',7r^ 



■ n^^{n) -H"^ = 0. 



(11) 



Suppose that m is a vector of length ek over GF(2''). We 
define the permutation vr of a vector performs only the column 



m 


Matrix 




inverse Fourier 






Transform 


c 




Multiplication 



















Fig. 1 . The block diagram of the encoding of QC codes in the transformed 
domain. 



permutations. Its reverse permutation is denoted by 7r~^. Then 
the vector c, defined as c = {c-^-'^}'^ ■il^^{n) = c-^-il^^{n) 
with 

c^^^ = m • G-^^'' (12) 

is named as transformed codeword. The vector c is a codeword 
of the QC code C from (fTTT i and (fT2] i. because 



= {m-G^^'^y' -n-^n) -11^ 
^ m- G^ -n-^in) -tl^ 
= 0. 



(13) 



Moreover, since 

|Q.F,7r|,r has full rank, the map- 

ping from m to c are one-to-one. As a result, the above 
equations can be viewed as an encoding in the transform 
domain. ETD is depicted in Fig. 1 and its steps are summarized 
as follows: 

Algorithm 1 The Encoding of QC Codes in the Transform 

Domain 

Input: 

The message m of ek symbols; 
The ek x en transformed generator matrix, G-^''^; 
Output: 

The binary transmitted codeword c of en symbols; 
Steps: 

1) The message m is encoded into the transformed code- 
word c-^ by the transformed generator matrix G-^''^ 

c-^ = {m • G-^-'^}'^"'. 

2) The transmitted codeword c is obtained by the inverse 
Galois Fourier transform from c-^'s. 



e-l 



(14) 



According the above encoding steps, it is straightforward 
to give the computational complexity of ETD. The first 
step of encoding involves e{n — k)k Galois additions and 
e{n—k){k — l) Galois multiplications. Since the computational 
complexity of the Fourier transform of length e implemented 
by the Galois fast Fourier transform (GFFT) [26] is less than 
e log2 e Galois operations, the computational complexity of the 
second step, which needs n GFFT, is less than ne log2 e Galois 
operations. Because log2 e is much smaller than [n ~ k) or k, 
the computational complexity of ETD is about 0{ek{n — k)). 

In terms of the memory consumption of ETD, it is clear that 
the biggest expenses on memory is for the transformed matrix 
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G-^ '^. Since it is a diagonal matrix, we only have to store e 
matrices, each of which requires the storage of (e{n — k)k) 
Galois symbols if it is systematic. In Table I, we compare 
the computational complexity and the memory consumption 
of ETD of QC codes with those of traditional encoding in 
terms of symbols. It is seen that the computational complexity 
of ETD is about TZ = e times lower than that of traditional 
encoding. 

Example 1. Consider a 6A-ary (4095, 2160) QC code with 
circulant size 63, i.e., n = 4095, k — 2160, n ~ k — 1953, 
e = 63 and r — 6. Thus, TZ = 63. In other words, the 
computational complexity of ETD is only 1.59% of that of 
traditional encoding. 

IV. Encoding of Binary QC Codes in the 
Transformed Domain 

In Section III, we derived ETD of QC codes over GF(2''). 
If the QC code is binary, both the message vector m and the 
codeword vector c are binary. However, there is no guarantees 
that the codeword c from ETD is binary. If there exist non- 
binary symbols in c, many more bits are needed to transmit 
it, resulting in code rate reduction. 

In this section, we propose a fast and simple mapping using 
subfield bases on the message vector m to make sure that the 
transmitted codeword c is also binary. In addition, to reduce 
the computational complexity of the first step of ETD, a post- 
processing step on the transformed codeword c-^''^, which is 
equivalent to pre-processing of the message vector m, will be 
presented. 

Recall the conjugacy constraint (|3]l for a binary vector. For 
the transformed codeword or c^'^, it means 



or 



'^ej + (2i)e 



-n(2i)e+J 



(15) 



(16) 



for < i < e and Q < j < n. For the sake of simplicity, 
we denote the transformed codeword and the codeword by n 
blocks of length e, c-^ — [c;^] and c = [cj], respectively. 



where c-f ,, — c-^. 



and Cj i — Cej+i- It is clear that Cj is the 



inverse Galois Fourier transform of cf. Similarly, we denote 



and 



by e blocks of length n, c 



[cf '"], where c. 



Then (flST l and (fTST i can be rewritten as 



(17) 



(18) 



respectively. 

We now show that if the message vector m is pre-processed 
by bases of subfields, then ([TtT i or (fTsT l is satisfied and hence 
the codeword c is binary. Again for the sake of simplicity, we 
denote the message vector m and the pre-processed message 
vector m by e blocks of size fc, m = [mi] and m = [rhi], 
respectively, where rrii^j = rriik+j and mi.j = rhik+j, i = 
0, 1, . . . , e — 1 and j = 0, 1, . . . , fc — 1. 

Recall the definition of the i-th conjugacy class = 
{Dt^, D°2, D°f' " } in ^ with size r]i. Suppose that a is 



a primitive element, then a basis /3i,0i A.ii • ■ • tPi.^-i spans 
the subfield GF(2''0 of GF(2'') whose element's 2''' power 
all equals to itself. Thus, if 7 is an element of GF(2'''), then 
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7 and "f — uiPi^i, where ui is in the ground field 



(=0 



GF(2). If rji = 1, then the subfield only has two elements, 
and 1. If rji = r, then the subfield is the field GF(2''). 

Using the bases f3i/s, we map message m to its pre- 
processed version lii via 



(19) 



First, the mapping from m to m is one-to-one since /3i,i's 
are linearly independent over GF(2). Second, m(2A'ti)e.j = 
"^(t") j' ^ detailed proof of this fact is given the appendix. 
The following theorem shows that the transformed codeword 



(20) 



encoded from the pre-processed message m satisfies the 
conjugacy constraint. 

Theorem 1. The transformed codeword t-'^ — {m.-G''^ ''^Y' 
encoded from the pre-processed message rh satisfies the conju- 



gacy constraint c^(2i)^ = (Cj^i)^ (US or c^:^_^ - = {cffY 
in ( I76l ). 

Proof: From the definition of c'^^'^, its {ni+j)-th symbol 
for i = 0, 1, . . . , e — 1 and < j < n is 

k-l 



(21) 



s=0 



From the definition of th (fT9] l and the conjugacy constraint 
on Dj in (|7ll, its (ri(2j)e + j)-th symbol is 



fc-1 



s=0 



s=0 



(22) 



Direct computation of (20) involves k Galois multiplications 
and fc — 1 Galois additions over GF(2'^), since both rh and 
are non-binary symbols over GF(2'"). Thus, we rewrite (21) as 



^{ii^ti) ,j ~ E "^(2''ti)^, 3-0(2(^4;)^, sj"j 

5=0 

fc-1 '- 



E E Pi,im(2H,)^,s 



s=0 \ 1=0 

J? — 1 /fc-1 

C E 

/=0 \s=0 



D 



2" 



E E -^H.us<s, 



(23) 



2'' 



1=0 ^ ' 

Equation ( |23] l shows that the mapping on message m fT% 
and the mapping on c-^ '^ ( |23] l result in the same c-^ '^, which 
satisfies the conjugacy constraint. 
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TABLE I 

Complexity of traditional encoding and transformed encoding 





Traditional Encoding 


Transformed Encoding 


Computational Complexity 


0{e'k{n~k)) 


0(efc((n-fc)) 


Memory Consumption 


ek{n — k) 


ek{n — k) 



However, the mapping on c'^''^ (l23T l involves much less 
Galois multiplications than the one on message m ( fT9] l. 
According to ( |23] |. computation of cf^'^ is carried out in 

two steps. In the first step, the summation of c^^i^ -^ ^ = 
fc-i 

Til ™(2'ti)c,s-D(2f ti) ,s,j is calculated, which only involves 

fc — 1 additions, since m is a binary vector; in the second 
step, cf^'" is calculated, which involves rji multiplications and 
rji — 1 additions. The other codeword symbol c^^^^^^ ■ in the 
conjugacy class of j can be simply calculated from (cf^'J)^^. 

Before we study the computational complexity of our pro- 
posed ETD algorithm, we summarize its encoding steps. For 
the sake of implementational simplicity, we assume that the 
permutation tt^^ operation is included in the initialization 
stage. 



Algorithm 2 The Encoding of Binary QC Codes in the 

Transform Domain 

Input: 

The message m of ek bits; 

The ek x en transformed generator matrix G-^''^; 
Output: 

The binary transmitted codeword c of en bits; 
Steps: 

1) The message m is encoded into the transformed code- 
word c-^ by the transformed generator matrix G-^ '^ 

c-^ = {m- G-^'^}'^"' 

2) The transformed codeword c-^ is mapped into the con- 
jugacy constraint satisfied codeword c-^, 

3) The binary transmitted codeword c is obtained by the 
inverse Galois Fourier transform from c-^'s, 

e-l 

Ci = ^ ^ C^ jOi\ 
3=0 



We present the computational complexity of the proposed 
ETD algorithm step by step. Since m is binary, the first step 
of encoding involves only ek{n — k) Galois additions. In the 
second step, each cjj. needs rji Galois multiplications and 
i]i — 1 Galois additions; the computational complexity of this 
step is thus about ne Galois multiplications and ne Galois 
additions. Since the GFFT over GF(2'') is about e log2 e Galois 
operations, so the complexity of the third step, which requires 



n times GFFT, is 0{ne\og2 e). 

In terms of the memory consumption, it is clear that most 
memory is spent on storing the transform matrix G"^ '^. Since 
it is a diagonal matrix, we only need to store e matrices, 
each of which requires to store {k{n — k)) Galois symbols 
if it is systematic. Moreover, considering the conjugacy con- 
straint O, we only need to store A representative matrices. 
Each symbol in the i-th representative matrix cost rji bits 
to store. Thus, the overall memory consumption of ETD is 

A 

T riikin — fc) = ek{n — fc), which is the same as that of 

'i=0 

traditional encoding. 

Considering that each Galois addition costs r bit operations, 
each Galois multiplication costs bit operations, and each 
Galois symbol costs r bits memory, we compare in Table II 
the computational complexity of ETD of binary QC codes with 
that of traditional encoding in terms of bits. 

Since r « log2 e. Step 1) costs most computation and 
dominates the computational complexity, which is about 
0{e{n — k)k logj e). Thus the complexity of the transformed 
encoding is about TZ = , ("~'')^'^ — = ^ — times lower 

^ (n — Kj/ce logo e ^og2 ^ 

than the complexity of the traditional encoding. 

Example 2. Consider a (4095, 2160) QC-LDPC with circulant 
size 63, i.e., n = 4095, k = 2160, n-k ^ 1953, e = 63 
and r = 6, we have, TZ ~ 10.05. In other words, the 
computational complexity of ETD is only 9.52% of that of 
traditional encoding. 

Example 3. Consider a (8176, 7154) QC-LDPC with circulant 
size 511, i.e., n 8176, k 7154, n - k ^ 1022, e = 
511 and r = 9, we have, TZ = 56.78. In other words, the 
computational complexity of ETD is only 1.77% of that of 
traditional encoding. 

V. Transformed Encoding of QC-LDPC Codes 

An LDPC code usually is defined by either its Tanner 
graph or its parity-check matrix. In this paper, we describe 
an {en,K) QC-LDPC code by its parity-check matrix, which 
can be simply denoted by a sparse (n ~ k) x n array of e x e 
circulants 

-A-o 1 ... Ao_„_i 
Ai.o Ai 1 ... Ai „_i 

H - . . . , 

A„_fe_i.o A„_/t_i^i ... A„_A;-l,r!-l 

(24) 

where K > ek. The authors of lfT2l proposed to compute the 
generator matrix G with quasi-cyclic structure from the parity- 
check matrix based on Gaussian elimination. Based on the ma- 
trix transformation (|6), we can have the transformed generator 
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Fig. 2. The block diagram of tiie encoding of QC codes in tiie transformed domain. 

TABLE II 

Complexity of traditional encoding and transformed encoding 





Traditional Encoding 


Transformed Encoding 


Computational Complexity 


0(e2fc(n-fc)) 


Step 1) 


ek{n — k)r 


Step 2) 


ne{r'' + r) 


Step 3) 


nr'^ e\og2 e 


Overall 


0{e{n — k)k log2 e) 


Memory Consumption 


ek(n — k) 


ek{n — k) 



matrix G-^ '^ from the generator matrix G. Alternatively, we 
can have G-^ '^ from the null space of the transformed parity- 
check matrix H'^ '^ [15|. In the sequel, we separately treat 
ETD of QC-LDPC codes with full rank H and rank deficient 
H. 



A. Encoding of QC-LDPC Codes with Full Rank Parity-Check 
Matrices 

Suppose that C is an {en,K) QC-LDPC code defined 
by its full rank parity-check matrix H = [A^j]. Sim- 
ilar to I©, we have the transformed parity-check matrix 
jjJ^.TT _ diagpo, Bi, . . . , Be_i]. Then we can construct 
the corresponding transformed generator matrix G"'^ '^ = 
diag[Do,Di,...,De_i] (with details in Us]) from H-^'P^ 
such that 



Pi = rankCBi) < n — k, i — 0,1, . . . , e — 1, such that 



= diag[JDoB'^, GiBJ, De-iBj_J = 0. 



To facilitate encoding, the diagonal matrices D^'s should 
satisfy the conjugacy constraint and be systematic. We thus 
propose to compute the systematic matrix from the matrix 
Bi, which is systemized from B^ by Gaussian elimination. In 

other words, if B^ = [P^ : I], then = [I : Pj] such that 



D,B7 = 0. Furthermore, D 



(D, 



\2f 



Suppose that H-^ '^ has full rank, i.e., K = ek. It is clear 

that B; = [Pi : I] for each i = 0, l,...,e — 1 has full rank 
and all B/s are of the same size. Then all the D^'s have the 
same size. Thus, ETD of QC-LDPC codes in this case is the 
same as before in Algorithms 1 and 2. 



B. Encoding of QC-LDPC Codes with Rank Deficient Parity- 
Check Matrices 

The parity-check matrices of many QC-LDPC codes, espe- 
cially algebraic QC-LDPC codes |5|, IHJ, [Ui, \M, lfT9l . 
II20I . II2TI . II27I . are rank deficient, i.e., K ^ ek. In some 
cases, more than half rows of their parity-check matrices are 
redundant. As a result, there exist rank deficient diagonal 
matrices in their transformed parity-check matrices. Conse- 
quently, their diagonal matrices B/s have different ranks, with 



B, 



P, : I 







(n—k) X n 

Thus, the diagonal matrices D^'s on the transformed gen- 
erator matrices have different sizes and ranks, with at = 
rank{T>i) — n — pi. The length of the message vector of the 

e-l 

QC-LDPC code is i^T = J2 '^i '^'^'^ '^^e transformed generator 

1=0 

matrix G^^'^ = [Dq, Di, . . . , De_i] can be represented by 



D, 



Thus, ETD of QC-LDPC codes with rank deficient parity- 
check matrices is similar to Algorithms 1 and 2 as well. The 
only difference lies in the first step because different D^'s 
have different ranks and sizes. Thus, ct^ symbols or bits are 
multiplied by in the first step in the non-binary or binary 



VI. Concluding Remarks 

In this paper, we have proposed a low-complexity encoding 
algorithm for QC codes in the transform domain. Its compu- 
tational complexity is much lower than traditional encoding, 
for both non-binary and binary QC codes. To further simplify 
ETD of binary QC codes, a post-processing step is devised to 
guarantee that the transmitted codeword is binary. In addition, 
it has been shown that ETD is applicable to both QC-LDPC 
codes with full rank parity-check matrices and QC-LDPC 
codes with rank deficient parity-check matrices. 

We have only considered cases with e = 2'' — 1 and q = 2'' 
or g = 2 in this paper. The derivations for other cases are 
similar. It is worth mentioning that, unlike traditional encod- 
ing, the transmitted codeword generated by our ETD algorithm 
is not systematic. Thus, after the codeword is corrected by 
the decoder of QC codes, Fourier transforms are required to 
recover the message vector The complexity of these Fourier 
transforms is much lower than that of traditional decoding 
algorithms for most QC codes. Moreover, non-systematic 
codewords by nature have better secrecy than systematic ones. 
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Appendix 



Theorem 2. Let rji be the size of i-th conjugacy class. The 
nonlinear mapping from a binary vector z to a vector z over 

GF(2''), 

Vt)e = (X] A^(2't)e)^^ (25) 
;=o 

is one-to-one correspondence, where {/?;} is a basis of the 
subfield GF(2'^^). Furthermore, z satisfies the conjugacy con- 
straint. 

Proof: First, we prove that it is a one-to-one correspon- 
dence mapping. It can be proved by contradiction. Suppose 
that there exist two distinct vectors z^^^ and z^^' which are 
mapped into the same vector z, i.e., 

("e A^^jj^^^fE (26) 

;=o ^ ' 1=0 ^ ' 

Since x"^ y"^ = (a; + y)"^ holds in the extension fields of 
GF(2), we have 

£ (4V - ^(2% J/5' = 0- 

Since z^^^ and z*^^^ are different, there exist nonzero co- 
efficients 2^^2i\j ~ • indicts that Pi's are linearly 
dependent, which contradicts the assumptions that Pi's are a 
basis. Similarly, it can be proved that any two different vectors 
z(^) and z'^) are mapped from two different vectors z^^^ and 
z(^^. As a result, the nonlinear mapping (25) is bijective. 

Then we prove that z satisfies the conjugacy constraint as 
follows. 

For < fi < r]i - 1, where r]i > 1, clearly, i(2M+it)^ = 
z^2(it) • For = ?7i — 1, we obtain that 

1=0 

Since /3; is over GF(2'''), we obtain /3f ' = Pi. Since z IS a 
binary vector, z^^H) ~ '^(2't)e- Thus 

Hence we can obtain that 

-2(2M+it)^ = -2(Vt)^> (27) 

for < < ?7i. 

Since (27) holds for all i, it follows that ^(2t)e = for 
< t < e. This is equivalent to the condition that z satisfies 
the conjugacy constraint. ■ 
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